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In this communication we consider the doping dependence of the strong antiferromagnetic spin 
fluctuations in the cuprate superconductors. We investigate the effect of an incommensurate mag- 
netic response, as recently observed in inelastic neutron scattering (INS) experiments on several 
YBa2Cu306+a; compounds, on the spin-lattice and spin-echo relaxation rates measured in nuclear 
magnetic resonance (NMR) experiments. We conclude that a consistent theoretical description of 
INS and NMR can be reached if one assumes spatially inhomogeneous but locally commensurate 
spin correlations and that NMR and INS experiments can be described within a single theoretical 
scenario. We discuss a simple scenario of spin and charge inhomogeneities which includes the main 
physical ingredients required for consistency with experiments. 
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I. INTRODUCTION 

Understanding the doping, frequency and tempera- 
ture dependence of the magnetic response in the high-T c 
cuprates is one of the most challenging problems for the 
high-T c community. Insight into the nature of the strong 
antiferromagnetic fluctuations very likely holds the key to 
the unusual normal state properties of the cuprates |],U , 
as seen in a variety of experimental techniques. Both 
nuclear magnetic resonance (NMR) and inelastic neu- 
tron scattering (INS) experiments arc important tools 
to probe spin excitations in cuprates While INS 

experiments provide insight into the momentum and fre- 
quency resolved imaginary part of the spin susceptibility, 
X(q, to), NMR experiments yield information on the mo- 
mentum averaged real and imaginary part of x(q, lj) in 
the zero frequency limit. However, in contrast to NMR 
experiments, INS measurements rely on the existence of 
large single crystals and often suffer from a rather limited 
experimental resolution. Therefore, these experimental 
techniques are complementary in the information pro- 
vided on spin excitations in the cuprates. 

One of the central questions in the cuprate supercon- 
ductors is whether INS and NMR experiments can be si- 
multaneously understood within a single theoretical sce- 
nario. Thus far it is not even been clear whether one can 
reach agreement between INS and NMR data as far as 
the order of magnitude of the spin susceptibility is con- 
cerned. This problem is caused, in part, by the fact that 
NMR is sensitive to all magnetic fluctuations, regard- 
less of whether they are related to a pronounced momen- 
tum dependence of the spin susceptibility. In contrast, 
INS typically probes only those parts of the susceptibil- 
ity which are strongly momentum dependent; the rest are 
typically attributed to the "background" of the signal, as 
might be caused by the scattering of neutrons on nuclei 
and phonons. 



The situation has been further complicated by re- 
cent INS experiments on YBa2Cu30g+a: which observe 
a crossover from a substantial incommensuration in the 
magnetic response at low frequencies, to a commensurate 
structure in \ at higher frequencies ^ g]. The question 
thus arises: does this incommensurate order originate 
from spatially inhomogeneous correlations between the 
doped holes |9| , which would preserve a locally commen- 
surate magnetic response, or does it reflect homogeneous 
incommensuration, most likely due to Fermi surface ef- 
fects. 

In this communication we argue that since NMR 
experiments probe the local spin environment around 
a nucleus, they are able to distinguish between a lo- 
cally commensurate and incommensurate magnetic re- 
sponse. Earlier theoretical studies of NMR experiments 
on YBa2Cu306+a; and YBa2Cu40s used a phenomeno- 
logical form of x PJ : 



x(q,w) 



l+£ 2 (q-Q) 2 -W^ s/ ' 



(1) 



where £ is the magnetic correlation length in units of the 
lattice constant, a , u> s { an energy scale characterizing 
the diffusive spin excitations, a an overall temperature 
independent constant, and Q is the position of the peak 
in momentum space which was assumed to be commen- 
surate, i.e. Q = (tt, it). Using this expression for x(Qi w ) 
a rather detailed quantitative understanding of various 
NMR data has been reached In particular, from 

the analysis of the longitudinal spin lattice relaxation 
rate of the 63 Cu nuclei, l/ 63 Ti, and the spin spin re- 
laxation time, 1/T 2 g, scann g l aws like ui s f oc £~ z with 
dynamical critical exponent, z, have been deduced. It 
was found that z«l between a lower crossover temper- 
ature, T*, and a higher one, T cr , whereas z ss 2 above 
T cr . In the temperature range where z = 1 scaling ap- 
plies, it follows that TiT/T 2 g is independent of temper- 
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ature; a conjecture which was experimentally verified by 
Curro et al. for YI^Ci^Os |14[. Below the pseudogap 
temperature, T», u> s i and £ decouple and a quasiparticle 
and spin pseudogap emerges. Moreover, it was recently 
shown that the temperature dependence of the incom- 
mensurate peak width, as determined in INS experiments 
on La2- 2; Sr 2 ;Cu04 ]15| , showed z — 1 scaling over a wide 
range of temperatures, in agreement with the NMR find- 
ings for YBa 2 Cu408- Because of its appearance in both 
NMR and INS results, we will therefore assume in the 
following analysis that z = 1 is the proper scaling be- 
havior for magnetically underdoped cuprates between T* 
and T cr . 

One of our central results is that the attempt to un- 
derstand the NMR data for YBa2Cu40s with homoge- 
neous incommensuration is inconsistent with z = 1 scal- 
ing. This implies that the local magnetic response, which 
is probed in NMR experiments, is commensurate. Should 
INS experiments show that this compound exhibits a 
(globally) incommensurate magnetic response similar to 
that seen in I^-rrSr^CuC^, this would strongly support 
a dynamic charge and spin inhomogeneity (stripe) origin 
of the incommensuration. In Sec. V we discuss a spin 
and charge inhomogeneity scenario which reconciles the 
global incommensuration with the local commensuration 
in x- 

In a second important result we carry out a quanti- 
tative comparison of the strength of the antiferromag- 
netic spin fluctuations, as measured by NMR and INS 
experiments. We find that agreement between the re- 
sults obtained from these quite different experimental 
techniques, which explore not only a different frequency 
range, but different wavevector regimes, can be obtained 
within a factor of 2. Given the large uncertainties in de- 
termining the absolute value of we believe that this 
result demonstrates that a consistent description of INS 
and NMR data can be achieved within the framework of 
Eq.(l). 

Although the form of x m Eq.(Q) was originally in- 
vented to understand the spin response at very low fre- 
quencies and above the superconducting transition tem- 
perature Jul], it has proved interesting to investigate to 
what extent one can understand INS data at higher fre- 
quencies within the same framework. The results by Aep- 
pli et al. support the picture of a unique incommensurate 
spin response from zero energy up to 15meV, the high- 
est energy used in Ref. |l5). Within the error bars of 
the experiment, only for momentum values away from 
the peak maximum do systematic deviations from Eq. (^) 
occur. Such deviations indicate the presence of lattice 
corrections to the continuum limit used in Eq. (jl|) . Lat- 
tice corrections are expected to be extremely important 
for the local, momentum averaged, susceptibility: 



(2) 



frequencies spin excitations away from the antiferromag- 
netic peak and energetically large compared to lo s { come 
into play and we therefore expect a poor description of 
x" oc (a;) in terms of Eq.(Q). This conclusion is indepen- 
dent of whether the peak position commensurate or in- 
commensurate peaks. 

Deviations from the universal continuum limit of x are 
expected to be also of relevance for the relaxation rates 
measured in NMR experiments since these are weighted 
momentum averages of x(q, to). Thus, it is worthwhile 
to study whether one can find indications for lattice cor- 
rections from an analysis of NMR and INS data. Such 
corrections to x are a l so of relevance for our understand- 
ing of the lifetime of so-called cold quasiparticlcs if one 
assumes that the lifetime of these quasiparticles is also 
dominated by scattering off spin fluctuation s ]16h 1S[] . 

Our paper is organized as follows. In Sec. [n[we give a 
brief overview of the theoretical framewor k in w hich we 
analyze the INS and NMR data. In Sec. [ill A 
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lyze NMR data on YBa2Cu40s for both commensurate 
and incommensurate magnetic res 
role of lattice corrections. In Sees. 



ponse and discuss the 
we an- 
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alyze NMR data on two YBa2Cu306+a; compounds. In 



Sec. IV we discuss the consistency between NMR and 



INS data, as well as the role of lattice corrections for 
the local susceptibility. In Sec. 0we propose a spin and 
charge inhomogeneity scenario as a possible way to rec- 
oncile NMR and INS data. Finally, in Sec. VI we draw 
our conclusions. 



II. THEORETICAL OVERVIEW 

We briefly discuss the theoretical framework in which 
we discuss NMR and INS experiments. In order to ana- 
lyze the low frequency NMR data, we use the dynamical 
spin susceptibility of Eq.(|l|), where for the commensu- 
rate case Q = (7r, tt). To allow for an incommensurate 
structure of x, we use Eq.(|l]) with Q = (1, 1 ± 5)ir and 
Q = (1 ± 5, l)ir and sum over all four peaks 11 1. The in- 
clusion of the correct form of lattice corrections is rather 
difficult since it requires a microscopic model which is be- 
yond the scope of this paper. In general we expect lattice 
corrections to appear in the form of an upper momentum 
cutoff, A. In the following we choose a soft cutoff proce- 
dure for <5q = q — Q in the denominator of Eq. ([!]) : 



<5q 2 



(3) 



since the phase space of momenta away from the peak 
maxima is considerable in two dimensions. For higher 



We also expect a weak momentum dependence of a and 
w s f once the continuum description breaks down; how- 
ever, to keep the number of tunable parameters small we 
ignore these effects. 

In NMR experiments, one measures the spin-lattice re- 
laxation rate 1/T\ X , with applied magnetic field in re- 
direction, and the spin-echo rate 1/T 2 g, which can be 
expressed in terms of the dynamical susceptibility as: 
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where a; = afe, c describes the direction of the external 
magnetic field. The 63 Cu form factors are given by 



63 F c (q) 
63 F e J f (l) 



A ab + 2B^ cos(q x ) + cos(q y ) 
A c + 2B\cos(q x ) + cos(q y )j 



Fab(q) 



r (>:i 



Ki f (q) + 63 F c (q) 



(6) 

(7) 
(8) 



where A a b, A c and B are the on-site and transferred hy- 
perfine coupling constants, respectively p| . 

It follows from Eq.(Q) that the spin excitations in the 
normal state are completely described by three parame- 
ters, a,£, and w s f. In order to extract these parameters 
from the experimental NMR data, we also need to ob- 
tain the three hyperfine coupling constants, A a b,A c , and 
B; hence we have six unknown parameters in the above 
equations, and require six equations to determine them. 
So far we have two, Eqs.(|4|) and (||). An additional con- 
straint arises from the temperature independence of the 
63 Cu Knight shift in a magnetic field parallel to the c 
axis m YBa 2 Cu 3 7 @ which yields 



A c + 4B « 



(9) 



A fourth constraint comes from the anisotropy of 
the 63 Cu spin-lattice relaxation rates |^(|, which for 
YBa 2 Cu307 was measured to be 
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3.7 ±0.1 



(10) 



A fifth constraint involving the hypcrfine coupling con- 
stants can be obtained by plotting the Knight shift 63 K a t, 
versus Xo(T) PJ, which yields 



AB + A n 
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A final constraint is obtained from the earlier conjec- 
ture that the antiferromagnetic correlation length at the 
crossover temperature T cr is approximately two lattice 
constants Jll| . Recent microscopic calculations have con- 
firmed this conjecture, showing that indeed £(T cr ) is of 
the order of a few lattice constants |2l[] . 

To the extent that experimental data for Ti and T 2 g 
are available up to T cr , one can fit the above set of 



six equations self-consistently to the data, and thus ob- 
tain not only the temperature independent parameters 
a,A ao ,A c , and B, but also the temperature dependence 
of £ and uj s f. It turns out, however, that this analysis 
is only possible for YBa2Cu40s, since this is the single 
compound for which data up to sufficiently high tempera- 
tures have been obtained |l4| ] . Our corresponding results 



will be presented in section III A . In order to extract 



the relevant parameters from NMR data on the related, 
widely studied, YBa2Cu306+a; compounds, and thus to 
make a comparison between NMR and INS experiments 
possible, we need to make two assumptions, which, at 
least partly, will be supported by the experimental data 
for YBa2Cu40g. First, we assume a relation between to s f 
and £, given by ||0 



uj sf = c£~ 



(12) 



where c is a temperature independent constant. For 
YBa2Cu40g, where we can independently extract Lu s f 
and £, we will show that there is indeed a crossover from 
z = 1 behavior at low temperature to z = 2 behavior at 
higher temperatures. The second assumption concerns 
the temperature dependence of the magnetic correlation 
length, £(T), which we take to be that obtained by one 
of us using a renormalization group (RG) approach [^2| : 
for z = 1 the temperature dependence of the magnetic 
correlation length is given by 



r 2 = Co 



bT 1 



(13) 



a result which is in good agreement with INS experiments 
on La 2 - a; Sr 2 ;Cu04 pl|. Having all the necessary theoret- 
ical tools in place, we now address the questions raised 
in the introduction. 



III. ANALYSIS OF NMR EXPERIMENTS 
A. YBa 2 Cu 4 8 

YBa2Cu408 is a benchmark system for our analy- 
sis, since Curro et al. |fl4|| have measured the relaxation 
times Ti and T2G over a wide temperature range from 
80K to 750K, and in particular between T*approx200 
K and T cr approx500 K. We are therefore able to ex- 
tract the relevant parameters which we discussed above 
from a self-consistent fit of Eqs.(Q)- ( ]TT| ) to the Ti and 
T2G data. In doing so, we assume that the constraints 
given by Eqs.(p|)-(pT|) for YBa 2 Cu307 are also valid for 
YBa2Cu408- Since we can determine £ and uj s f inde- 
pendently, we are able to study the effect of both an in- 
commensurate magnetic response, and of non-universal 
lattice corrections, on the scaling law uj s f — c/£ z . 

We first reconsider the case of a commensurate struc- 
ture of x, and set A = 2y / 7r, i.e. we choose a momentum 
cut-off which corresponds to the linear size of the BZ. We 
present the temperature dependence of both w s j and £ , 
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FIG. 1. The temperature dependence of u s f and in 
YBa2Cu40g for the commensurate case with A = 

which results from a self-consistent fit to the experimen- 
tal data, in Fig. Between T, = 200 K and T cr = 500 
K, uj s f and £ _1 clearly scale linearly with temperature. 
In order to study the extent to which u s f and £ obey the 
above scaling law, we plot in Fig. g, ln(w s /) as a function 
of ln(£) (lower curve). The dynamical scaling range is 
of course too limited to prove the existence of a scaling 
relation. However, assuming u s { oc £ -z , we hnd z w 1 
below 500 K, which we identify with T cr and z ~ 2 above 
T cr . These results are consistent with an earlier analysis 
of the scaling behavior |lj,[l3| and the microscopic sce- 
nario of Ref . . To study the effect of stronger lattice 
corrections, we decreased A to A = 2y / 7r/4 (upper curve 
in Fig. H which for clarity is offset). Our conclusions con- 
cerning the scaling relations remain unchanged. We can 
thus conclude that the dynamical scaling behavior of spin 
excitations is robust against even sizeable lattice correc- 
tions, as long as the structure of x" 1S commensurate. 

We next examine the effect of incommensuration on 
the dynamical scaling. It was earlier argued that the 
magnetic response in YBa2Cu40s should be very sim- 
ilar to YBa2Cu306.8 Using the observed doping 
dependence of the incommensuration in YBa2Cu306+x, 
we estimate the incommensuration in YBa2Cu40s to be 
S = 0.23. Setting A = we plot in Fig. || ln(w s /) as 
a function of ln(£) (lower curve) and find, following the 
same argumentation as above, that the dynamical scaling 
exponent below T cr sa 500 K is now z ps 0.75. Increas- 
ing the strength of the lattice corrections by decreasing 
A we find that z is increased. However, in order to ob- 
tain again z fa 1 (the upper curve in Fig. [|), which one 
would expect from the INS data on La2- a; Sr 2 ;Cu04, we 
need a very small momentum cut-off, A ss 2y / 7r/15. Such 
a small cut-off of order 0(l/£) is inconsistent with the 
continuum theory which is the basis for a dynamical scal- 
ing approach. We therefore conclude that a spatially ho- 
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FIG. 2. ln(u> s f) as a function of In^) for the commensu- 
rate case. Upper data set: A = 27r 1//2 /4; Tower data set: 
A = 2n 1/2 . Solid line: z=l; dashed line: z=2. 

mogeneous incommensurate magnetic response is in con- 
tradiction with the available NMR data. The physical 
origin of this sensitivity of the magnetic response of a 
homogeneously incommensurate system, compared to a 
locally commensurate one, results from the fact that the 
former effectively decreases the spatial extent of the spin- 
spin correlations, thus increasing the role of large values 
of Sq ~ A. It should also be noted that our arguments 
using lattice corrections to discriminate between these 
two scenarios works only for intermediate values of the 
correlation length. For very large £ the system should be 
insensitive to the cut-off procedure regardless of whether 
it is commensurate or incommensurate. 

In what follows we assume that the low frequency mag- 
netic response measured in an NMR experiment is indeed 
locally commensurate which implies that z = 1 scaling 
prevails between T* and T cr . In Sec. |y| we present a pos- 
sible theoretical scenario to resolve this apparent contra- 
diction between the incommensuration seen in INS exper- 
iments and our conclusions based on the available NMR 
data. Finally, using the commensurate form of Eq. ([!]), 
and a momentum cut-off A = 2t/it we find the following 
parameters from the solution of the above self-consistent 
equations: a = WeV^ 1 , A a b = 20.2fcOe/^s, A c — 
-182.8 kOe/nB, and B = 45.7 kOej[i B for YBa 2 Cu 4 8 . 

B. YBa 2 Cu 3 7 

Since it was earlier argued that for the optimally doped 
YBa 2 Cu 3 07, T cr ta 125 K, only slightly above T c , there 
is no significant temperature dependence of the relax- 
ation rates between T c and T cr . We will therefore per- 
form our analysis of the NMR data only at T cr . Assum- 
ing that £(T cr ) = 2, and setting A = 2-^7?, we utilize 
T 2G data by Itoh et al. || and Stern et al. ||, as 
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FIG. 3. ln(uj 3 f) as a function of Zn(£) for the incommen- 
surate case. Lower data set: A = Upper data set: 
A = 2tt 1/2 /15. Solid line: z=l; dashed line: z=0.75. 



well as Tic data by Hammel et al. Q25J to find the fol- 
lowing hyperfine coupling constants A a b = 27.8 kOe/fj,B, 
A c = -175.2 kOe/(i B , and B = 43.9 kOe/^ B - Note that 
these values agree well with the hyperfine coupling con- 
stants we extracted for YBa2Cu40s. With 1/T 2 g = 10 4 
s _1 and T lc T w 0.15 Ks at T cr , we obtain a = 18.5 eV _1 , 
u s / 19.8 meV and c « 39.6 meV. It turns out that the 
above parameter set is rather robust against changes in 
A. When A is decreased from A = 2y/n to A = 20F/4, A c 
and B decrease by about 1.5%, A a b increases by about 
8.5%, a increases by 11%, and oj s f decreases by about 
15%. A change of the momentum cut-off by a factor of 
4, which decreases the area of integration in the BZ by 
a factor of 16, thus leads only to moderate changes in 
the parameter set. However, as we show in Sec. IV, the 



corresponding changes in the local susceptibility at high 
frequencies are much more dramatic. 



C. YBa 2 Cu 3 06 63 

In what follows we show that, though only a limited 
set of NMR data on YBa2Cu3 06.63 is available, we never- 
theless reach the same conclusions regarding the scaling 
behavior as in YBa2Cu40s- In particular, since no NMR 
data on YBa 2 Cu306.63 are available above T = 300 K, 
and since we do not know the exact value of T cr , we can- 
not extract the parameter set from a self-consistent fit 
to the NMR data as we did in the case of YBa2Cu40s- 
However, it is in this doping range that INS data are 
available and therefore a consistency check between the 
parameter sets extracted from NMR and INS measure- 
ments might be the most promising. It turns out that 
though we are not able to perform a fully self-consistent 
fit, we can still extract the parameter set if we make sev- 
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TABLE I. Parameter sets for two different values of T c , 
and A. 



eral additional assumptions which are supported by our 
results on YBa 2 Cu40s. First, we assume that z = 1 
scaling is present between T* and T cr and that we can 
describe the temperature dependence of £ and the re- 
lation between ui s f and £ by Eqs.(|l2|) and (|l3|), re- 
spectively. Second, we assume that the hyperfine cou- 
pling constants are only weakly doping dependent, and 
that to good approximation we can use the same con- 
stants we extracted from the analysis of YBa 2 Cu3C>7 for 
YBa2Cu3C>6.63- Third, we need an estimate for T cr , 
which is unknown for YBa 2 Cu3 06.63- However, since 
T cr w 500 K for YBa2Cu40s, and since T cr increases 
with decreasing doping, we assume T cr w 550 — 650 K for 
YBa2Cu3C>6.63- To demonstrate the effect of the uncer- 
tainty in the latter assumption, we calculate he parame- 
ters a, uj sf {T) and £(T) for T cr = 650 K and T cr = 550 K 
as well as for two different values of the momentum cut- 
off A. The resulting parameter sets based on experimen- 
tal data by Takigawa et al. p6| a re shown in Table |. Note 
that the parameters in Table III are rather robust against 
large variations in T cr and/or A and only differ at the 
most by about 30%. In Fig. | we compare our theoreti- 
cal results for the relaxation rates, 1/Ti and l/Tgc, with 
the experimental data by Takigawa et al. |2q| . We find 
that the temperature dependence of the relaxation rates 
calculated for each of the four different parameter sets in 
Table |, is practically indistinguishable; we therefore only 
plot the results for the parameter set with T cr = 650 K 
and A = 2y / 7r. This result is consistent with the ro- 
bustness of the scaling behavior in the commensurate 
case against non-universal lattice corrections as found for 
YBa 2 Cu 4 08. We find a consistent description of the ex- 
perimental data for YBa2Cu3 06.63 using the assumption 
of a z — 1 scaling relationship and local commensuration 
in the temperature regime « 230 K < T < 300 K, and, 
as was the case for YBa 2 Cu40s, assuming commensurate 
behavior leads to results which are practically indepen- 
dent of non- universal lattice corrections. 
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FIG. 4. Comparison of our theoretical results for the tem- 
perature dependence of 1/XiT (upper figure), and I/T2G 
(lower figure) in YBa2Cu30e.63 with the experimental data 
by Takigawa et al. [M. 



IV. INELASTIC NEUTRON SCATTERING 

In this section we address the following two ques- 
tions: (i) What is the effect of lattice corrections on the 
local susceptibility, and (ii) can one describe the avail- 
able INS data with the parameter set extracted in the 
previous section ? 

As we discussed in the introduction we expect lattice 
corrections to lead to substantial corrections of x" oc at 
frequencies larger than u> s f. Above this energy scale we 
also expect that the typical energy, A sw , of propagat- 
ing spin waves, which, at low frequencies are completely 
overdamped by particle hole excitations, comes into play, 
leading to a modified form of the spin susceptibility 



1 + £ 2 (q - Q) 2 - iu/u s} - (w/A 8W ) 2 



(14) 



For uj < Ag W /cj s f no sign of a propagating peak exists 
due to the diffusive character of the spin excitations. For 
oj > Ag W /u s f and q = Q the consequence of a propa- 
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FIG. 5. The local susceptibility Xi'od 1 ^) as a function of 
frequency for two different values of A. 



gating mode is a pronounced pole in the excitation spec- 
trum if A sw < uj s f and a soft upper cut off in energy if 
A sw > uj s [. The form of x(q, oJ) in Eq.(|l4|) can be shown 
to describe the propagating spin mode in YBa2Cu30g.5 
above T c 

Using the parameter sets (1) and (2) in Table B, ex- 
tracted in the previous section from NMR experiments, 
we plot in Fig. g the local susceptibility of YBa2Cu30g.63 
for two different values of A. We clearly see that upon 
decreasing the momentum cut-off, the maximum in Xi oc 
moves towards lower frequencies. Experimentally, how- 
ever, the intensity at higher frequencies drops much faster 
than in Fig. |L even for A = 2v / rF/4 plfl . In order to 
explain this discrepancy, we note that at higher frequen- 
cies, the dominant contribution to x'Lc comes horn re- 
gions in momentum space which are far away from the 
peak position at (tt, tt). In these regions, x"(*l; ^-0 is on ly 
weakly momentum dependent in which case its contribu- 
tion might be easily attributed to the experimental back- 
ground. In other words, we believe that the origin of the 
discrepancy lies in an underestimate of the experimental 
intensity at higher frequencies due to the problems one 
has in resolving it from the background (29[ . 

In Fig. |^ we plot the calculated INS intensity, i.e., 
X"(Q,w) as a function of frequency at Q = (tt, tt). Since 
the calculated intensity is practically the same for all four 
parameter sets in Table we only present x"(Q; w ) calcu- 
lated with the second parameter set in Table | (solid line). 
In the same figure, we also include the experimental data 
by Fong et al. for YBa 2 Cu306.7 J|(|. For the tempera- 
ture range of interest, this material is the closest match 
for YBa2Cu3 06.63- Calculating the overall INS intensity 
with the parameters extracted from NMR experiments in 
the previous section, we find that we underestimate the 
experimental INS intensity by about 56 %, or a factor of 
2.3 (the dashed line shows the calculated INS intensity, 
multiplied by an overall factor of 2.3). Given the uncer- 
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taintics in the experimental determination of the absolute 
scale of x"j we believe that the above result represents 
reasonable agreement between the INS and NMR data 
and thus demonstrates consistency in the description of 
spin excitations based on these two experimental tech- 
niques. 



V. SPIN AND CHARGE INHOMOGENEITIES 



We saw in Sec. Ill A that under the assumption of 
z = 1 scaling, supported by recent INS experiments [ fl5"| , 
a locally incommensurate magnetic response is inconsis- 
tent with the available experimental NMR data. We 
therefore concluded that NMR data support a locally 
commensurate magnetic structure. The question thus 
arises of how one can understand a locally commensu- 
rate, but globally incommensurate magnetic response as 
seen by INS? It has been suggested that charge stripes 
separated by an average distance = 2tt/S are the origin 
of the incommensurate magnetic response seen in INS ex- 
periments 1 3 1 1 . For a large part of the sample this would 
leave the locally commensurate structure intact. Except 
for the Nd-doped 214 compounds B^] and probably for 
the particular doping value x = ^~ in I^-ajSr^CuC^, 
these stripes are believed to be dynamic rather than 
static in nature. 

Considering a spatial charge variation: 



p(r) = Po + Sp(r) , 



(15) 



where po is independent of r. The question arises whether 
the system is in a regime with p Sp(r), as suggested in 
Ref. fnf, or rather characterized by rather small charge 



variations, i.e. 5p(r) < pq. The particular appeal of the 
former assumption of strong charge modulations is an 
immediate explanation of the doping dependence of the 
position of the incommensurate peak, 5 = 2x. Stripes 
with one hole every second lattice site in a hole free an- 
tiferromagnetic environment move closer together upon 
doping, leading to the above x-dependence of 8. Also, 
the exceptional behavior of various magnetic and trans- 
port properties for the doping value x = | can be easily 
understood in terms of a stable commensurate arrange- 
ment of such stripes and the underlying lattice. Another 
argument in favor of this scenario are the recent results 
by Vojta and Sachdev |3^] who investigated the mean 
field behavior of a system with competing magnetic and 
long range Coulomb interactions and found areas with 
a strong charge modulation separated by regions with 
barely disturbed antiferromagnetic correlations. 

Nevertheless, there are several conceptional problems 
with such a strong charge modulation, which lead us 
to present some arguments favoring a more moderate 
charge variation, 5p(r) < po, under circumstances that 
the transverse mobility of the charge carriers with respect 
to the averaged position of the stripes is large. First, 
uncorrelated spatial stripes of width ro, which fluctuate 
transversely over a distance d, will give rise to a broaden- 
ing d/lg of the magnetic peaks observed by INS. The ex- 
perimental INS results, however, provide strong support 
for a scenario in which the width of the magnetic peaks 
is determined by Jig] . Thus, in order to observe sep- 
arated incommensurate peaks, we need at least £ d to 
establish a well defined antiphase domain wall, implying 
d <ti Iq. In other words, assuming uncorrelated stripes, 
the width d over which stripes fluctuate must be very 
small to account for the existing INS data. Such " stiff' 
but uncorrelated stripes seem to be in contradiction to 
the notion of stripes as a dynamical entity. Therefore, we 
arrive at the conclusion that spatial stripe fluctuations, 
if they exist, must be strongly correlated. 

Second, at the characteristic temperatures where 
stripes with po ^ $p( r ) appear, strong modifications 
of the resistivity and other transport properties are ex- 
pected to occur. None or only moderate changes of this 
kind have been observed and no indication for a dimen- 
sional crossover from quasi one-dimensional to quasi two- 
dimensional dynamics seems to be present. 

Third, for low frequencies, the spin excitations in 
doped cuprates are overdamped rather than propagating 
spin waves. The suppression of the spin damping upon 
opening of a quasiparticle gap in the superconducting as 
well as in the pseudogap state implies that the dominant 
source of the spin damping, 7, are particle-hole excita- 
tions. Assuming weak charge modulations, we then find 
Imx _1 (w) oc 70?, in agreement with the results of INS 



experiments 15 1. 

On the other hand, rigid stripes separated by one di- 
mensional charge carriers, which are effectively bosonic 
in character, lead to strong deviations from the above 
frequency dependence of Imx _1 (tj) p4|,|35||, in disagree- 
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merit with INS experiments. Thus, Imx _1 (o;) does not 
arise from spatially varying stripes even though they cer- 
tainly affect the spin degrees of freedom and can cause 
the decay of these excitations. 

Finally, NMR and other magnetic measurements on 
various cuprates which provide strong support for spatial 
inhomogeneities |36|-f4l|] suggest mostly spatially varying 
spin degrees of freedom but not necessarily a strong mod- 
ulation of the charge background. 

The above points suggest that the inhomogeneities ob- 
served in the high-T c cuprates might be predominantly 
magnetic in origin with only moderate modulations of the 
charge density. This is not unrealistic since in strongly 
correlated systems small variations in the charge den- 
sity can bring about substantial changes in the magnetic 
properties of the system. This scenario of a predomi- 
nantly magnetic character of the inhomogeneities, lead- 
ing to magnetic stripes and clusters, is schematically pre- 
sented in Fig. ^. Note, that the magnetic stripes, in 
which the electronic charge density is lower than that in 
the magnetic clusters, represent magnetic domain walls 
with a phase slip of 7r in the ordering of the spins. Theo- 
retically, we believe that this more homogeneous charge 
arrangement is a result of quantum fluctuations beyond 
the mean field level investigated in Ref. (33| . 

The formation of magnetic clusters and stripes requires 
competing interactions on different length scales. Be- 
sides the local Coulomb interaction, we need to identify 
a longer length scale interaction. It was argued earlier 
that in the cuprate superconductors an increase of the 
magnetic correlation length leads to strong vertex cor- 
rections, which in turn give rise to a gradient coupling 
between the fermionic and spin degrees of freedom p2[ . 
For a system with an ordered ground state, this takes the 
form of a dipole-dipole coupling which has been shown 
to cause inhomogeneities and domain formation in two 
dimensional magnetic systems. The length scale of this 
gradient coupling should roughly speaking be the mag- 
netic correlation length, £. We thus need a minimum 
value of £ to (a) generate the gradient coupling, and to 
(b) observe the incommensuration. The temperature at 



which incommensuration should appear is roughly set by 
£(T) w Iq = 4, where the last equality arises from the mo- 
mentum position of the incommensurate peaks. Since, as 
we argued earlier, £(T cr ) 2, T cr thus presents an up- 
per bound for the formation of magnetic clusters. This 
anomalous coupling is of course enhanced once magnetic 
clusters are formed; the creation of magnetic inhomo- 
geneities is a self-consistent process. The magnetic in- 
homogeneity scenario presented here is admittedly very 
qualitative and further microscopic investigations will be 
required to verify or disprove it. 

VI. CONCLUSIONS 

In this communication, we investigated whether it is 
possible to understand theoretically INS and NMR data 
within a single theoretical framework. Based on the ob- 
servation of an incommensurate structure of the mag- 
netic response by INS experiments |^,^5|, we considered 
whether this reflects a locally or globally incommensu- 
rate ordering. On analyzing NMR data on YBa2Cu40s 
we found with the condition of z — 1 scaling that a ho- 
mogeneous incommensuration is, within the framework 
given by Eq.(|l]), inconsistent with the available experi- 
mental data. We thus concluded that the local magnetic 
structure as probed by NMR is likely commensurate. 

We then investigated the effect of lattice corrections 
on both the parameter sets extracted from NMR data 
and on the local susceptibility measured in INS exper- 
iments. We found that while the resulting corrections 
to the NMR parameters are only weak, the most pro- 
nounced effect of such corrections appears in the local 
susceptibility determined in INS measurements at fre- 
quencies above lo s $. Specifically, we found that decreas- 
ing the momentum cut-off leads to a suppression of x'i oc 
at higher frequencies, thus improving the agreement with 
the experimental data. We argued that the remaining 
discrepancies can be explained by an experimental un- 
derestimate of the INS intensity at higher frequencies. 
Furthermore, we quantitatively compared INS and NMR 
data and found agreement within a factor of 2. Given the 
large uncertainties in resolving x" from the large back- 
ground and in determining the absolute intensity of \" , 
we believe that this result demonstrates that a consistent 
description of INS and NMR data can be obtained using 
the expression for \ given in Eq. (|l]) . 

Finally, we discussed a spin and charge inhomogeneity 
scenario to reconcile the local commensurations, as seen 
in NMR experiments, and the incommensurate peaks, 
seen by INS. This scenario, even though similar in spirit 
to earlier charge stripe pictures, is based on a moder- 
ate to weak modulation of the charge density, causing 
a pronounced inhomogeneity in the magnetic properties 
and leading to magnetically coupled clusters separated 
by weakly correlated stripes. 

We would like to thank A.V. Chubukov, P. Dai, B. 



8 



Keimer, T. Mason, H. Mook, R. Stern, CP. Slichter and 
B. Stojkovic for valuable discussions. This work has been 
supported in part by the Science and Technology Center 
for Superconductivity through NSF-grant DMR91-20000 
(D.K.M.), and by DOE at Los Alamos (D.P.). 



[1] D. Pines, Z. Phys. B103, 129 (1997); Proc. of the NATO 
ASI on The Gap Symmetry and Fluctuations in High T c 
Superconductors, J. Bok and G. Deutscher, eds., Plenum 
Pub. (1998). 

D. J. Scalapino, Phys. Rep. 250, 329 (1995). 
C. P. Slichter, in Strongly Correlated Electron Sys- 
tems, ed. K. S. Bedell et al. (Addison-Wesley, Reading, 
MA, 1994). 

T. E. Mason, bond-mat/9812287j 

Ph. Bourges, |cond-mat/9901333| ; Ph. Bourges, Y. Sidis, 
H. F. Fong, B. Keimer, L. P. Re gnault, J. Bossy, A S. 
Ivanov, D. L. Milius, I. A. Aksay, sond-mat/9902067 



J.M. Tranquada et al, Phys. Rev. B 46, 55 61 (1992) 



P. Dai et al,. preprint, cond-mat 9712311 



H. Mook et 



al., preprint, cond-mat 9712326. 

P. Dai et al, Phys. Rev. Lett. 80, 1738 (1998). 

Due to the absence of a corresponding elastic Bragg peak 

in neutron and X-ray scattering, stripes, if present at all 

in YBa2Cu306+:r, have to be dynamical in character. 

A. Millis, H. Monien, and D. Pines, Phys. Rev. B 42, 

1671 (1990). 

Y. Zha, V. Barzykin and D. Pines, Phys. Rev. B 54, 7561 
(1996). 

A. Sokol and D. Pines, Phys. Rev. Lett. 71, 2813 (1993). 
V. Barzykin and D. Pines, Phys. Rev. B 52, 13585 (1995). 
N. J. Curro, T. Imai, CP. Slichter, and B. Dabrowski, 
Phys. Rev. B 56, 877 (1997). 

G. Aeppli, T. E. Mason, S. M. Hayden, H. A. Mook, J. 
Kulda, Science 278, 1432 (1997). 

R. Hlubina and T.M. Rice, Phys. Rev. B 51, 9253 (1995); 
ibid 52, 13043 (1995). 

B. P. Stojkovic and D. Pines, Phys. Rev. Lett. 76, 811 
(1996); Phys. Rev. B 55, 8576 (1997). 

L. B. Ioffe and A. J. Millis, Phys. Rev. B 58, 11631 
(1998). 

B. Stojkovic, private communication 
S. E. Barrett, et al. Phys. Rev. Lett. 66, 108 (1991). 
J. Schmalian, D. Pines, and B. Stojkovic, Phys. Rev. 
Lett. 80, 3839 (1998); i bid to appear in Ph ys. Rev. B. 



:ond-mat/9810041 



J. Schmalian, preprint, 
Y. Roh et al, J. Phys. Soc. Jpn 63m 1455 (1994). 
R. Stern, M. Mali, J. Roos, and D. Brinkmann, Phys. 
Rev. B 51, 15478 (1995). 

P.C Hammel et al, Phys. Rev. Lett 63, 1992 (1989). 
M. Takigawa et al, Phys. Rev. B 43, 247 (1991). 
D.K. Morr and D. Pines, Phys. Rev. Lett. 81, 1086 
(1998). 

P. Dai, et al, Science 284, 1344 (1999). 
P. Dai, private communication. 



[30] 
[31] 



[32] 
[33] 
[34] 
[35] 
[36] 

[37] 
[38] 



[39] 

[40] 

[41] 
[42] 



H.F. Fong et al. Phys. Rev. B 54, 6708 (1996). 
See, for example, J. Zaanen and O. Gunnarson, Phys. 
Rev B 40, 7391 (1991); V. J. Emery and S. A. Kivelson, 
Physica C 263, 44 (1996); O. Zachar, S. A. Kivelson and 
V. J. Emery, Phys. Rev. B 57, 1422 (1995). 

J. M. Tranquada, cond-mat-/9802043. 

M. Vojta and S. Sachdev, preprint, cond-mat/9906104 . 

A. H. Castro Neto, Z. Phys. B103, 185 (1997). 
J. Schmalian {unpublished). 

J. H. Cho, F. C. Chou, and D. C. Johnston, Phys. Rev. 
Lett. 70, 222 (1993). 

F. Borsa et al, Phys. Rev. B 52, 7334 (1995). 

B. J. Suh, P. C. Hammel, Y. Yoshinary, J. D. Thom- 
son, J. L. Sarrao, and Z. Fisk, Phys. Rev. Lett. 81, 2791 
(1998 ).;P. C. Hammel , B. J. Suh, J. L. Sarrao, and Z. 



Fisk, |cond-mat/980909e . 

A. W. Hunt, P. M. Singer, K. R. Thurber, and T. Imai, 
Phys. Rev. Lett. 81, 5209 (1998). 

M.-H. Julien, F. Borsa, P. Caretta, M. Horvatic, C. 



:ond-mat/9903005 



Berthier, and C. T. Lin, 
J. Haase et al. preprint. 

JR. Schrieffer, J. Low Temp. Phys. 99, 397 (1995); A.V. 
Chubukov and D.K. Morr, Phys. Rep. 288, 355 (1997); J. 
Schmalian, D. Pines, and B. Stojkovic, Phys. Rev. Lett. 
80, 3839 (1998). 







